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We investigate the linear-chain structures in highly excited states of ^"^C using a generalized 
molecular orbital model, by which we incorporate an asymmetric configuration of three a clusters in 
the linear-chain states. By applying this model to the system, we study the ^°Be-|-a correlation 
in the linear-chain state of ^''C. To clarify the origin of the ^"Be-l-a correlation in the linear- 
chain state, we analyze linear 3a and "ia+n systems in a similar way. We find that a linear 3a 
system prefers the asymmetric 2aH-a configuration, whose origin is the many body correlation 
incorporated by the parity projection. This configuration causes an asymmetric mean field for two 
valence neutrons, which induces the concentration of valence neutron wave functions around the 
correlating 2a. A linear-chain structure of '^^C is also discussed. 
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p I. INTRODUCTION 

The linear-chain structure of three a clusters in C isotopes has been a long standing problem, for example, as 
I . discussed in Refs. [ujlL Historically, the 0^ state in ^^C was suggested to be a candidate for the linear-chain structure 
] of three a clusters [jQl • However, according to 3a cluster model calculations in the 1970s, it was pointed out that 
' ^ '■ the state is not a linear-chain state, but a weak coupling state, where three a clusters are weakly interacting 
■ with each other [^-Q. In recent years, the state attracts an interest as the a condensed state in analogy to 
,_H , Bosc Einstein condensation in nuclear matter [lol - [l2 |. Instead of the state, the state was suggested to have a 
^ ■ geometric configuration of three a clusters [isj . In the recent works on ^^C with antisymmetrized molecular dynamics 
2 \ (AMD) and fermionic molecular dynamics (FMD) calculations [l4| - [T6j . the Oj^ state was suggested to have an open 
I triangle structure of three a clusters and was assigned to the experimental 0"*" state at 10.3 MeV. However no linear- 

chain structure of three a clusters nor linear-chain band was found in the ^^C system, and therefore, searching for a 
O^l , linear-chain state in C isotopes is still an attractive subject. 

Linear-chain structures in C isotopes have been discussed in earlier works with dia+Xn cluster models 0-11, [l3| 
and the Skyrme mean- field model In the study by Itagaki et al. [1], they assumed a symmetric 3a cluster 
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■ configuration with valence neutrons in molecular orbitals. It was argued that the linear-chain structure of ^''C may 
! be unstable because there is no local minimum in the energy curve against the bending mode. 

■ However, to confirm the stability of an excited state, it is essential to take into account the orthogonality to lower 
I states. In contrast to the work of Ref. in our previous study on ^''C with an AMD method [l^, we found 

the linear-chain structure and its rotational band consisting of the 0'^ , 2^, and 4^ states in ^^C. In Ref. [l^, we 
adopted a method of f3-j constraint AMD in combination with the generator coordinate method (GCM) called the 
/3-7 constraint AMD + GCM [l^. Although the 3a core is not a priori assumed in the framework, various 3a-|-2n 
structures are obtained in the energy variation. As a result of energy variation and GCM calculations, the linear-chain 
states appear and show ^''Be-l-a correlation. Here, the ^''Be-l-a correlation means an asymmetric 2a-|-a configuration 

■ with two valence neutrons distributing around the correlating 2a. We note that three a clusters have one geometric 
configuration, and therefore ^''Be-l-a (2a-|-a) correlation does not have a weak coupling feature, but a strong coupling 
feature. In other words, this ^"Be is not an isolated one, but a polarized configuration of the middle a cluster and 
valence neutrons in the linear-chain structure. 

For excited states of "'^^C, there exist many experimental works suggesting the ^"^Be-l-a cluster structure in the 
observed excited states above the ^"Be-l-a threshold [13, [20l - l23| . We expect that these states might be possible 
candidates for the linear-chain states. Then we consider that the ^^Be+a correlation is one of the key features of the 
linear-chain structure. 

Our aim is to confirm the ^"Be-l-a correlation in the linear-chain states of ^^C. As discussed in Ref. [l^, our previous 
work indicates that the stability of the linear-chain structure in ^*C is owed to the valence neutron effects and the 
orthogonality to lower states. Since the bending mode is suppressed by the orthogonality in a ^*C system, a 3a 
cluster core with a straight line configuration is considered to be rather stable against the bending mode. Therefore, 
in the present study, we consider only the linear configuration and omit the bending mode to find the origin of the 
^"Be-l-a correlation. For ^^C, this ansatz is probably valid. In such a linear system, the ^"Be-l-a correlation can be 
described by the asymmetry of three a clusters and valence neutron motion. Here the first question is whether three 
a clusters favor an asymmetric 2a-|-a configuration which were not incorporated in the work by Itagaki et al. 
The second question is whether the two neutrons localize and form the ^'^Be cluster or if they are moving around all 
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three a clusters as described in the molecular orbital picture [171 ] . Unfortunately, in our method used in the previous 
work [l^ , detailed behavior of valence neutrons in the molecular orbitals is not sufficiently described because a single 
particle wave function is approximated by a Gaussian wave packet. 

In the present paper, to investigate the linear-chain structures of "'^^C we extend the molecular orbital model by 
explicitly incorporating asymmetric configurations of the linear-chain structure. We solve the valence neutron motion 
in the linear 3a system by using a linear combination of p orbits around a clusters as done in the molecular orbital 
model. By analyzing the calculated results, we will confirm the ^^Bc-t-o; correlation in the linear-chain structure of 
^^C and clarify its origin. To see the effect of three a clusters and that of valence neutron motion, we compare the 
results of ^^C with those of 3a and 3a-|-n. Here, we comment that these 3a and 3a-|-n do not necessarily correspond 
to excited states of realistic ^^C and -^'^C systems. Although the Q'^ state of ^^C is suggested to have an open triangle 
structure, it is different from the straight-line chain structure assumed in the present analysis. Even if a linear-chain 
structure does not exist in realistic ^^C and "'^'^C, analysis of 3a and 3a-|-n systems with the linear-chain configuration 
is useful to understand the origin of ^°Be-|-a correlation suggested in the realistic linear-chain states in ^^C. We also 
calculate ^^C and discuss systematics of correlations in the linear-chain structure of C isotopes. 

This article is organized as follows. In SecHIl we explain the framework of the generalized molecular orbital model. 
The calculated results are shown in Sec. IIIII We discuss the motion of valence neutrons and correspondence to the 
previous work in Sec. IIVI Finally, in Sec. |Vl we give a summary. 



II. FRAMEWORK OF THE GENERALIZED MOLECULAR ORBITAL MODEL 



In this section, we explain the framework of the generalized molecular orbital model. In this model, we express 
the valence neutron motion in the linear-chain structures by a linear combination of configurations written by the p 
orbits centered at a clusters. 



A. Wave function of the generalized molecular orbital model 



In the present model, the wave function of a carbon isotope with a mass number A, ^C, is described by a super- 
position of basis wave functions which consist of three a clusters with a linear-chain configuration and X ^ A — \2 
valence neutrons with various configurations: 



(1) 



where A is the antisymmetrizing operator and c represents a set of the configuration for valence neutrons, c = 
{ci, • • • , cx}- The subscripts L, M, and R mean left, middle, and right positions for three a clusters, respectively. We 
label the left, middle, and right a clusters aL, cum, and aR. The coefficients /(c) are determined by diagonalization 
of Hamiltonian and norm matrices as explained later. 



The wave function of at is given by the (Os) configuration whose center is shifted at a certain position Rt 
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where t = L, M, R. A linear-chain structure is expressed by the positions Rt on the z-axis and it is characterized by 
inter-cluster distances, di and d2- 



di 
d2 



\Rm-Rl\ 
\Rr-Rm\ 



as seen in Fig. [TJ Without loss of generality, we can limit the distances d2 > di. 

For spatial wave functions of valence neutrons, we define the p orbits centered at one of three a clusters, 

{P±)t = -^i{px)t ±i{Py)t) , 

{Px)t = (l>Rt+Se^ - (t>Rt-Se^, 
{Py)t = 4>Rt+Sey - (pRt-Sey, 



(4) 
(5) 



(6) 

(7) 
(8) 
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FIG. 1: The schematic figure for a hnear-chain structure of three a clusters. The distances between a clusters L and M, and 
M and R are set to d\ and d2, respectively. 



where 5 is a small value. The {p±)t represents the p orbits around at, and it has the z-component of orbital angular 
momentum = ±1. In general, the fc-th valence neutron wave function ^p^J.{ck) in Eq. [T]is given by a product of {p±)t 
and a spin wave function, 

<Pfc(cfc) = (p±)tXnt or {p±)tXni- (9) 

(p+)tXrit and (p-)tXni are the eigenstates of the z-component of total angular momentum with j'^ = SI = +3/2 and 
VL ~ —3/2, respectively, which are energetically favored because of the spin-orbit force, while {p+)tXni and {p-)tXnt 
are the unfavored fl = ±1/2 states. In the present model, we omit the p^ orbits because the orbits constructed from 
the pz orbits have higher nodes than those from p^ and Px orbits; therefore the Pz orbits are unfavored. 

In the present calculation, we assume the first and second valence neutrons occupy the favored orbits and the third 
and fourth valence neutrons occupy the unfavored ones. As a result, the valence neutron wave functions can be written 
with the labels n = ±3/2, ±1/2 and t = L, M, R as, 

ip^h) = {p+)t,Xn^ = [+3/2]t„ (10) 
¥^2(^2) - {P-h^Xni ^ [-3/2]t„ (11) 
¥'3(i3) = (p+)t3Xnt=[+l/2]t3, (12) 

V2iU)^{p-)uXni^[-mt,, (13) 

and a set of the configuration for valence neutrons, c, can be explicitly written as c = {ti, • • • ,tx} with ti = L, M, R. 
It means that the total wave function is a superposed wave function given by the configuration mixing for valence 
neutrons with respect to ti = L, M, R indicating positions of the p orbit centers, 

'i'-c^ E fih,--- ,tx)A{ipl,^^,^^,ip'iih),--- ,^3citx)}. (14) 

The coefRcient f(ti,--- ,tx) is determined by the diagonalization to minimize the energy of the state. for a 

given 3a configuration parametrized by di and ^2. 

In this model, wc introduce an asymmetric correlation of a spatial configuration of a clusters by the position of the 
middle a cluster. When the middle a cluster shifts to one side considerably, d2 ~ di ^ fm, a linear three a system 
has the correlating 2a and the isolated a configuration. This improvement is the most different point from symmetric 
molecular orbital models. Moreover, we incorporate a multincutron correlation by the configuration mixing of valence 
neutron orbits. 



B. Parity and angular momentum projections 

We project the wave functions onto parity and angular momentum eigenstates by using the parity projection 
operator and the angular- momentum projection operator P'lu^. The parity projection operator is defined as 

where P is the parity operator. The angular-momentum projection operator -P//^ is defined as 

_ 2J + 1 /■ J, 



PuK ^ I dnDi^K{n)R{n), (i6) 
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where Q = {a, (3,^) are the Euler angles, D'ljj^{^) is the Wigner's D function, and R{fl) = e~'"'^^e~''''^"e~'^'^"' is the 
rotation operator. 

We will show energy curves of intrinsic, parity projected, and parity and angular-momentum projected wave func- 
tions. 

C. Interaction and parameters 

For the effective two-body interactions, we use the Volkov No. 2 interaction [2^ as the central force and the 
spin-orbit term of the GSRS interaction [25| as the LS force. We take the same interaction parameters as those in 
Rcfs. i.e., the Majorana exchange parameter M — 0.6 {W = 0.4), the Bartlett exchange parameter B = 0.125, 

and the Heisenberg exchange parameter H = 0.125 in the central force, and ui = —1600 MeV and U2 = 1600 MeV in 



the LS force. These parameters are the same as those adopted in the studies for ^Be l26l . and ^°Be [27[, except for a 
small modification in the strength of the spin-orbit force to fit the 0^ energy of ^^C |la |. 

For the width parameter of single-particle Gaussian wave packets in Eq. ([3]), we used the value ly = 0.235 fm~ 
which is also the same as those in the studies for C isotopes 0, EB, E^, HI] . 



D. One Slater determinant wave function 



To discuss the correlation between valence neutrons in a linear-chain structure, which may be incorporated by the 
diagonalized wave function with the configuration mixing in Eq. 1141 we also introduce one Slater determinant wave 
function having no explicit two-body correlation in principle, and compare the results. 

In the framework of one Slater determinant wave function, each single particle wave functions for a valence neutron 
is given by a linear combination of [57] t orbits with respect to t = L, M, R instead of the configuration mixing in Eg. 1141 

$o-si-t" = ^{^«,^^,<^g,^"(f]i),--- ,7r-(17x)}. (17) 

where 

Tr-in)= Cn,tMt- (18) 

t=L,M,R 

ill = +3/2, ^2 = —3/2, ila = +1/2, and SI4 = —1/2 are chosen for the first, second, third and fourth valence neutrons 
as well as Eg . [Til and C^^t arc the coefficients for the orbits {^]t- n~{il) means the character of the orbit: the tt orbit, 
negative parity, and z component of angular momentum f2. In the limit of Cq.l = Cq^m = C'n.R = 1 and di ~ c?2, 
these TT~ (Q) orbits coincide with the tt orbits of a usual molecular orbital model. 

In the present study, the coefficients C^.t are determined by a variational calculation. In a variational calculation, 
we apply the following conditions for the coefficients, 

C+3/2,t = C^3/2.t, (19) 

C+i/2,t = C'-i/2,t- (20) 

The model space spanned by this one Slater determinant wave function is a subspace of the generalized molecular 
orbital model. 

We also introduce three kinds of test wave functions with fixed coefficients which are expected from naive pictures. 
We name three orbits 'MO', 'PO', and 'AO'. In 'MO', the coefficients are chosen as 

Cn,L = Cq^m = Cn,R = 1, (21) 

which indicates the molecular orbital limit where valence neutrons spread over the whole of three a clusters. 'AO' 
corresponds to the atomic orbital limit with valence neutrons localizing around the middle a cluster: 

Cn,M — 1, (22) 
Cn,L = Co,R = 0. (23) 

'PO' means the partial orbit in which valence neutrons distribute around the correlating 2a as 

Cn,L = Co,M = 1, (24) 
Cn.R = 0. (25) 



FIG. 2: (Color online) The energy surfaces of positive-parity states in ^"'C, positive-parity states in ^^C, and negative-parity 
states in ^''C. The a-a distances, di and d2, are taken as described in Fig. [l] 

This distribution describes a naively expected Be correlation. However, please note that realistic Be correlation in the 
linear-chain state of ^^C is different from this 'PO' configuration as shown later. For example, the 7r~ (-1-3/2) orbits 
for 'MO', 'PO', and 'AO' arc written as follows: 



Neutrons in 'MO' orbits gain the kinetic energy because valence neutrons distribute in a wide region along the z-axis, 
while those in 'AO' orbits gain maximum potential energy, as shown later. 



We applied the generalized molecular orbital model to C isotopes: "'^^C, ^^C, "'^'^C, and ^^C. In this section, we show 
the calculated results. 

The energy surfaces of positive-parity states in ^"^C, positive-parity states in ^^C, and negative-parity states in ^'^C 
are shown as functions of ^2 — d\ and d\ ~\- d2 in Fig. [2l di + d2 means the total length of the linear chain, while d2 — di 
is the difference of the intervals between a clusters and it indicates asymmetry of the 3a configuration (see Fig. [ij . 
Namely, linear three a clusters have symmetric structure for d2 — di ~ Q fm, and the asymmetric structure develops 
as the d2 — di value increases. At the minimum point of the energy surface of ^^C, the linear-chain length c?i -I- d2 
and the difference d2 — di are 5.2 fm and 1.5 fm, respectively. The existence of the energy minimum at the finite 
d2 — di value indicates ^^Be+a correlation in ^"'C, though this minimum is shallow. At the minimum points of the 
energy surfaces of ^^C and ^'^C, the differences d2 — di are 1.6 fm and 1.5 fm, respectively, which are similar to the ^^C 
case. However, the linear-chain lengths di + d2 are 6.5 fm and 5.8 fm, respectively, which are different from the ^^C 



'MO' : 7r-(+3/2) 
'PO' : 7r"(+3/2) 
'AO' : 7r-(+3/2) 




[+3/2]l + [+3/2]m + [+3/2]r 



[+3/2]m. 



(26) 
(27) 
(28) 



III. RESULTS 
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case. More valence neutrons shorten the hnear-chain length but do not change the difference of the intervals between 
a clusters. Although the locations of energy minima are different depending on the number of valence neutrons, all 
the energy surfaces indicate Be+a correlation in the linear-chain structures of C isotopes. Here, we stress again that 
the aim of calculating the ^^C and ^^C systems with linear 3a and Sa+n structures is to clarify the origin of the 
^"Be+a correlation in ^^C. In the following, we fixed the length di + d2 to be 5.2 fm and show the energy curves of 
the linear-chain structure as a function of the difference d2 — di because the energy surface of ^^C is rather soft over 
the wide range of the difference d2 — di but is steep in the linear-chain length di + d2 around the di + d2 = 5.2 fm 
line. 

In Fig. |3l we show the energy curves of ^"'C, ^^C, and ^^C with the linear-chain structure. In each panel, the 
energy curves for the intrinsic, parity projected, and angular-momentum projected wave functions of the generalized 
molecular orbital model are shown. To discuss asymmetry of the linear-chain structure, we shift the energy curves 
for the intrinsic and parity projected wave functions of ^''C, ^^C, and ^^C by —6.9 MeV, —6.1 MeV, and —6.2 MeV, 
respectively, to adjust the energies to the angular-momentum projected energies at d2 — di = fm. 

At first, we explain the results for ^^C shown in the top panel of Fig. [3l The intrinsic energy curve has the energy 
minimum at the symmetric configuration, d2 — di = fm, while the positive-parity and the 0"*" energy curves have 
energy minima at asymmetric configurations, d2~di > fm. This difference indicates that the many body correlation 
incorporated by the parity projection gain the energy with the asymmetric configuration of three a clusters resulting 
in the 2a+a correlation. The kinetic energy gain by the parity projection, which spreads the center of mass motion of 
a clusters, becomes larger than the interaction energy loss in the asymmetric configuration at the energy minimum. 
The depth of the 0+ minimum, 0.73 MeV, is deeper than that of the positive-parity minimum, 0.56 MeV. Moreover 
the position of the 0"*" minimum, d2 — di = 1.7 fm, is more asymmetric than that of the positive-parity minimum, 
d2 — di = 1.5 fm. That is, many body correlation incorporated by the angular- momentum projection enhances the 
the asymmetric correlation in the present calculation. 

Next, we discuss the results for ^^C (the middle panel of Fig. |3]) to extract the contribution of three a clusters 
without valence neutrons to the asymmetric correlation. The behaviors of the energy curves are qualitatively similar 
to that of ^^C. The positive-parity and 0"*" energy curves have minima at asymmetric configurations, while the intrinsic 
energy curve has the minimum at d2 — di = fm of the symmetric configuration. The depth of the O'^ minimum, 1.02 
MeV, is deeper than that of the positive-parity minimum, 0.67 MeV. Again, the degree of the asymmetry of the O'^ 
energy curve {d2 — di — 2.4 fm) is larger than that of the positive-parity energy curve (d2 — di = 2.0 fm) . Interestingly, 
even if valence neutrons arc absent, a 3a system prefers an asymmetric configuration resulting in the 2a correlation 
after the parity projection. That is, linear three a clusters have the *Be-|-a correlation, and ^'^Be+a correlation in the 
linear chain states of originates in this nature of 3a systems. Compared with the parity projected and O"*" energy 
curves for ^^C with those for -'^''^C, it is found that energy curves in the large d2 — di region are less steep, and the 
energy minimum points are deeper and more asymmetric than those for ^^C. The former is easily understood from 
the viewpoint of Be-a inter-cluster potentials. That is, the potential energy between ^''Be and a clusters is deeper 
than that between ®Be and a clusters because of the additional two neutrons. 

From the results for ^^C shown in the bottom panel of Fig. [31 we can see the effect of a valence neutron in the linear- 
chain structure. The behaviors of the energy curves are qualitatively analogous to those of ^'^C and ^^C. Although 
the intrinsic energy curve has the minimum at the symmetric configuration, the negative-parity and the 3/2~ energy 
curves have minima at asymmetric configurations. The depth of the 3/2~ minimum, 0.79 MeV, is deeper than that 
of the negative-parity minimum, 0.57 MeV. The degree of the asymmetry of the 3/2^ energy surface {d2 — di = 2.0 
fm) is larger than that of the negative-parity energy surface (^2 — rfi = 1.7 fm). 

In the systematics of the minima of the angular momentum projected energy curves for ^^C, ^'^C, and ^^C, it is found 
that the asymmetry c?2 — di is the largest in the linear 3a structure without valence neutrons. As a valence neutron is 
added to the linear 3a system one by one, the depth of minima becomes shallow and the asymmetry c?2 — c?i reduces 
slightly. Therefore valence neutrons seem to somehow reduce the asymmetry of the linear 3a structure. However, as 
we show in the following, valence neutrons gather around the middle of linear three a clusters and play an important 
role to enhance the Be-f a correlations. 

We discuss the motion of valence neutrons and their effect in the linear 3a structure by comparing the results of the 
generalized molecular orbital model with those calculated by using simpler wave functions, 'MO', 'PO', 'AO', and one 
Slater determinant defined in Sec. IIIDI In Fig.|4l we show the angular-momentum projected energy curves of ^'^C and 
^•^C as a function of the difference ^2 — di of the a-cluster intervals. The top panel shows the ^^0(0''') energy curves 
for the 'MO'; TO', 'AO', and one Slater determinant wave functions compared with the generalized molecular orbital 
ones, and the bottom panel shows the ^■^C(3/2~) energy curves. It should be noted that one Slater determinant wave 
function is equivalent to the generalized molecular orbital wave functions in -'^'^C, because two models are coincident 
with each other in the case of a single valence neutron. This is not necessarily the case in ^''^C because the generalized 
molecular orbital wave functions contain two neutron correlations due to the configuration mixing, while one Slater 
determinant does not have two neutron correlations. 
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FIG. 3: Energy curves of C, C, and C with the linear-chain structure whose length di + 6,2 is fixed to be 5.2 fm, as 
a function of the difference of the a-cluster intervals, 6,2 ~ d\. Each panel shows the energy curves for the intrinsic, parity 
projected, and angular-momentum projected wave functions of the generalized molecular orbital model. The intrinsic and 
parity projected energy curves are shifted to agree with the angular-momentum projected energy curves at di — d2 = fm. 



In spite of the difference of the model space, the obtained one Slater determinant wave functions for ^ C give the 
same 0"*" energy curve as that of the generalized molecular orbital model in all regions of 6,2 — di. Actually, the 
obtained ^^C wave functions of these two models are found to be equivalent to each other with about 99% overlap. 
This means that the two-body correlation between valence neutrons is not essential in the linear chain structure of 
^^C but the intrinsic wave functions can be described by a simple product of single-particle wave functions as written 
by one Slater determinant in Eq. [TT] Next, we describe the energy curves obtained by using three kinds of test wave 
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FIG. 4: Energy curves of C and C with the hnear-chain structure whose length is fixed to be 5.2 fm. Energies are plotted 
as a function of the difference of the «-cluster intervals, d2 — di. The top and bottom panels show the ^'*C(0^) and ^'^C{3/2^) 
energy curves for the 'MO', 'PO', 'AO', and one Slater determinant wave functions compared with the curve for the generalized 
molecular orbital model wave function. 



functions, 'MO', 'PO', and 'AO'. Although the energy curves are higher than those of the generalized molecular orbital 
model and one Slater determinant, all of these wave functions show soft energy curves against the asymmetry d2 — di 
similarly to the case of the generalized molecular orbital model. The energy minima of 'PO' and 'AO' reproduce an 
asymmetric configuration. In particular, the degree of the asymmetry of 'AO', d2 — di = 1.9 fm, is in good agreement 
with that of the generalized molecular orbital model, though the depth of the 'AO' energy minimum is larger. 

We discuss valence neutron wave functions in the one Slater determinant wave function of ^^C, which almost 
coincides with that obtained by the generalized molecular orbital model. As mentioned before, the coefficients C3/2,t 
[t = L, M,R) in Eq. [T7]are optimized for each ^2 — di. At the energy minimum, d2 ^ di = 1.7 fm, the ratio of the 
coefficients is 

C3/2,L : C^3/2,M : C3/2.R = 1.0 : 4.3 : 0.6. (29) 

The largest coefficient C3/2,m uidicates that valence neutrons gather around the middle a cluster(Q;M)- Since the au 
is shifted to one side, the concentration of the amplitude at the q;m means that valence neutron distribution leans 
toward the narrow side of two of three a clusters. That is, the asymmetric configuration of linear three a clusters 
and valence neutrons following the middle a cluster make ^"Be+a correlation, which is consistent with the result in 
Ref. [l9| . For more quantitative discussion of the asymmetry, the reduction of the positive parity component in the 
intrinsic wave function can be a good measure. At the d2 — di = 1.7 fm corresponding to the energy minimum of 
linear-chain states in ^'^C, the magnitude of the positive parity component in the linear 3a system without valence 
neutrons is 89%, while that with two valence neutrons for the whole ^^C system is found to be 82%. The reduction of 
the positive parity component, i.e., the enhancement of the asymmetry by valence neutrons owes to the concentration 
of the neutron distribution at the middle a cluster. Incidentally, 82% for ^^C at the energy minimum, c?2 — di = 1.7 
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FIG. 5: The positive-parity energy curves of C using 'MO', 'PO', 'AO', one Slater determinant, and the generalized molecular 
orbital model with the linear-chain structure whose length is 5.2 fm. 

fm, is very close to the magnitudes of the positive parity component, 81%, in the linear-chain structure of ^^C at the 
energy minimum, d2 — di = 2.4 fm. In the next section, we show that this localization of valence neutrons can be 
understood by a mean field potential from three a clusters with the linear-chain structure. 

In ^^C, the situation is quite similar to ^"'C. Although the energy curves of 'MO', 'PO', and 'AO' are higher than 
that of the generalized molecular orbital model, they give energy minima at asymmetric configurations. As well as 
^^C, the degree of the asymmetry of 'AO' at the energy minimum agrees well with that of the generalized molecular 
orbital model. Here, we also analyze the distribution of the valence neutron in ^'^C. At the energy minimum, the ratio 
of the coefficients 6*3/2,4 (.7 = L, M, R) is 



As well as ^"'C, the valence neutron gathers around the middle a cluster and the ^Be+a correlation is found in ^^C. 

Surprisingly, the ratio of coefficients and the behaviors of energy curves in ^''C are almost the same as those in ^'^C. 
This is consistent with the fact that the linear can be described by one Slater determinant wave function and 
valence neutrons have no explicit dineutron correlation. In this analysis, we find that valence neutrons in are 
moving rather independently in the mean field constructed from the linear three a clusters. In other words, a valence 
neutron in ^^C is not affected by another valence neutron. 

Finally, we show the energy curves of ^^C to study the systematics of correlations in linear chain structures of C 
isotopes. For ^^C, we do not calculate the 0+ energy curve because of the computational cost. However, we have 
already seen in the results of i2,13,i4q w^qx the 0"*" energy curves are essentially the same as those of the positive- 
parity energy curves. Therefore, in this paper, we discuss the positive-parity energy curves of ^^C. In Fig. [5l we 
show the positive-parity energy curves of one Slater determinant and the generalized molecular orbital model wave 
functions. The behaviors of the energy curves are different from those of other C isotopes. The most striking feature 
is that the positive-parity energy curve has no energy minimum at an asymmetric configuration. Although it is soft 
around d2 — di =0 fm, an asymmetric correlation seems to disappear in ^^C. Moreover, in ^^C, the energy curve 
of the generalized molecular orbital model is lower by 1 MeV than that of one Slater determinant wave function in 
the d2 — di < \ iva region. This means the configuration mixing in the generalized molecular orbital model gains 
some correlation energy due to multi- neutron correlation in the linear-chain structure of ^®C. Comparing the wave 
function at the symmetric configuration d2 — di ^ fm between the one Slater determinant model and the generalized 
molecular orbital model, the overlap with each other is found to be 95%, which is still large, but smaller than the 
99% overlap for the case of ^''C. 

The ratios of the coefficients Cn,t (fi = 3/2, 1/2, t = L, M, R) in the one Slater determinant wave function are 



As same as ^^C and ^^C, the valence neutrons gather around the middle of the linear a clusters. However, in ^^C, 
the 3a core does not show the ^^Be-|-a correlation but it is symmetric at the energy minimum. Furthermore, the 
valence neutrons also have the symmetric configuration because they distribute around the middle a cluster, which 



C3/2X ■ C3/2M ■ C'3/2,R = 1.0 : 4.4 : 0.8. 



(30) 



6*3/2, L ■ 6*3/2,14 
Cl/2,L ■ C*i/2,M 



C*3/2,R = 1.0 : 4.9 : 0.5, 
C*i/23 = 1.0 : 5.0 : 0.7. 



(31) 
(32) 
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z[fm] 

FIG. 6: (Color online) The single-folding potential of three a clusters with the linear-chain structure. The crosses show the 
positions of three a clusters. 

exists just at the center of linear-chain structure. This result indicates a linear-chain structure in ^^C has a symmetric 
configuration. 

The question to be answered in the future is whether or not the linear-chain structure is stable in excited states of 
^^C. Firstly, the formation of three a clusters should be checked by such a framework free from model assumptions 
of clusters as the AMD method pol - is^ . Secondly, because the orthogonality to low-lying states is essential for the 
stabilization of linear-chain states as found in ^''C [l^, it is necessary to investigate the structure of ^^C covering 
low-lying states to high-lying states systematically. 



IV. DISCUSSION 



In this section, we give further discussion of the behavior of valence neutrons in the linear-chain structure. As 
already mentioned, valence neutrons tend to distribute around the middle a cluster. It seems in contrast to a simple 
molecular orbital picture, where valence neutrons are considered to distribute widely surrounding all a clusters. The 
tendency of the concentration around the middle a is understood by a mean-field potential formed by linear three a 
clusters as follows. To see the effect of the mean-field potential for a valence neutron from linear three a clusters we 
calculate the single-folding potential using the Volkov No. 2 force, 

U{r) = j drV(r')^'™''"'(r' - r'), (33) 

of three a clusters with the linear-chain structure. Here p{r') is the density distribution of linear three a clusters. 
In Fig. [6l we show the contour plot of the folding potential on the xz plane. Valence neutrons feel this potential in 
linear three a clusters. As clearly seen, the folding potential is deepest in the narrow side of two of three a clusters 
and becomes shallow around the isolated a cluster; therefore, valence neutrons tend to gather around the narrow 
side to gain potential energy. Around the narrow side, valence neutrons feel attractions from the two a clusters and 
gain larger potential energy, while around the isolated a cluster, valence neutrons feel an attraction only from one a 
cluster and gain less potential energy. 

Here, we note that valence neutrons do not distribute the deepest region of the potential because of Pauli blocking 
with a clusters. In fact, neutron wave functions in p orbits, {p±)t (t = L,M,R), vanish on the x = y = axis. In 
the present calculations, the root mean square distance \/{p^ from the x = y ^ Q axis for a p orbit with the width 
parameter v ~ 0.235 fm~^ is 



\/- = 2.06 fm. (34) 



Roughly speaking, valence neutrons in p orbits feel the potential around the y{p^ = 2.06 fm surface of the linear 
three a structure. The folding potential along the x = 2.0 fm line is shown in Fig. [T] The behavior of this energy 
curve is consistent with that seen in Fig. [B] That is, the folding potential is deepest in the narrow side of two of three 
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FIG. 7: The folding potential of three a clusters with the linear-chain structure for x = 2.0 fm. The crosses show the positions 
of three a clusters 




FIG. 8: energy surfaces of C with the linear-chain structure against the difference between the distances of q clusters. The 
energy surface for the mixing configurations of 'MO', 'PO', and 'AO' are illustrated as well as that of the generalized molecular 
orbital model. 



a clusters and becomes shallow around the isolated a cluster. The mean field constructed from linear three a clusters 
induces the concentration of valence neutron wave functions around the middle a cluster. 

As shown above, the mean-field potential induces the concentration of valence neutron wave functions around the 
middle a cluster, while the localization of the wave functions loses kinetic energy. As a result of the balance of potential 
energy gain and kinetic energy loss, valence neutron wave functions arc optimized. Here, we describe features of three 
kinds of test wave functions, 'MO', 'PO', and 'AO', and discuss valence neutron motion from a viewpoint of kinetic 
and potential energies. In 'MO', valence neutrons spread over the whole of three a clusters. 'PO' describes a naively 
expected Be correlation with valence neutrons distributing around two of the three a clusters. 'AO' means valence 
neutrons around the middle a cluster. 'MO' gains larger kinetic energy because of the spread of valence neutron wave 
functions, while it gains less potential energy because of no concentration around the middle a cluster. On the other 
hand, 'AO' gains larger potential energy and looses kinetic energy because of the localization. The failure to reproduce 
the results of the generalized molecular orbital model by the 'PO' model indicates that the Be correlation in linear 
three a clusters is different from the naively expected Be correlation which is incorporated in the 'PO' wave function. 
In the obtained results of the generalized molecular orbital model, both the kinetic energy gain and potential energy 
gain play important roles. To demonstrate this feature, we show the 0+ energy curves of ^^C calculated by mixing two 
of the three 'MO', 'PO', and 'AO' configurations in comparison with the results of the generalized molecular orbital 
model in Fig. [51 The point is that the 'MO+AO' curve is very close to that of the generalized molecular orbital model. 
This indicates that the linear-chain state in -'^"'C can be interpreted as an intermediate state of 'MO' and 'AO' to gain 
both of the kinetic energy and potential energy. On the other hand, other sets of mixing, 'MO+PO' and 'PO+AO' 
fail to reproduce the results of the generalized molecular orbital model and 'MO-I-AO' mixing. 
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Next, we discuss the relation between the present results and our previous results with AMD [l^- In the present 
results, the minimum point of the 0+ energy surface (^2 ~ di = 1.7) has ^°Be+a correlation and 80.4 MeV binding 
energy. In the AMD results, the main component of the 0^ linear-chain state also has ^^Be+o; correlation and 80.6 
MeV binding energy. The binding energies are almost the same and the valence neutron behaviors are qualitatively 
similar. Moreover, the overlap between these wave functions is significant as 81%. Therefore, we can conclude that 
the present results correspond to the AMD results. This good correspondence between the present and AMD results 
indicates that the previous AMD calculation properly describes the valence neutron configurations in the linear-chain 
structure of the 0^ in ^^C. In Ref. [l9|, we expected that the linear-chain states (rotational band consisting of the 
0^, 2q, and 4^ states) might be the candidates for the states observed recently in ^^Bc-f-o; breakup reactions. The 
^*^Be-|-a correlation confirmed in the present study supports again this expectation. 

We comment on the quantum fluctuation of the a-a distances. The full solution does not correspond to well defined 
distances but to a superposition of the wave functions which have various distances because the energy curves are very 
soft with respect to asymmetry and quantum fluctuation may be large. However, because of the existence of minimum 
at the asymmetric point, components of asymmetric configuration are expected to largely mix in the superposed wave 
functions. It means the importance of Be-|-a correlation. 

Finally, we comment on the negative parity states of the linear-chain structures in ^"^C. Because of the asymmetry 
of the intrinsic state from the ^"Be-l-a correlation, one can expect negative parity states as a parity partner of the 



positive-parity linear-chain states |17| . However, in the present calculations, energies of negative parity states are 



much higher, more than 10 MeV, than the positive-parity linear-chain states; therefore, the negative-parity linear- 
chain states can not appear near the positive-parity linear-chain states. It is consistent with the previous study on 
^^C with an AMD method [l^ which shows a negative-parity linear-chain state does not appear. 



V. SUMMARY AND OUTLOOK 



We investigated the linear-chain structures of ^ C using a generalized molecular orbital model. In this model, we 
incorporated asymmetric configurations of linear three a clusters to study the ^"Be-l-a correlation in ^''C. We also 
investigated Be+a correlations in linear 3a and 3a+n systems to see the effect of three a clusters and that of valence 
neutron motion. The linear-chain structure of ^^C is also studied. 

It was found that the linear-chain structure of ^^C shows the ^°Be-fa correlation, which is characterized by the 
asymmetric 2Q!-t-a configuration with two valence neutrons distributing around the correlating 2a. Compared with 
the results for the linear 3a and 3a+n systems, we find that a linear 3a system prefers the asymmetric 2a-|-a 
configuration, whose origin is the many body correlation incorporated by the parity projection. This configuration 
causes an asymmetric mean field for two valence neutrons, which induces the concentration of valence neutron wave 
functions around the correlating 2a. These are the origins of the ^"Be-l-a correlation in the linear-chain structure of 
"'^^C. In analysis of motion of valence neutrons, we find that both the kinetic energy gain and potential energy gain 
play important roles in the linear three a clusters with valence neutrons. A symmetric configuration of linear three 
a clusters is indicated in ^^C. We also find valence neutrons in a linear-chain structure of ^^C have multi- neutron 
correlation, while those in ^^C have no explicit two-neutron correlation. 

As a future problem, it is important to clarify whether or not a linear-chain structure with a symmetric configuration 
appears in ^^C. For this aim, the AMD -I- GCM method is a promising framework for systematic study of ^^C from 
low-lying states to high-lying states. 
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